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Structure functions of the nucleon and nuclei in neutrino reactions 



S. Kumano ^ * 

'^Department of Physics, Saga University, Saga, 840-8502, Japan 

Determination of parton distribution functions in nuclei is important for calculating nuclear corrections in 
oscillation experiments, from which detailed information should be extracted on neutrino properties. First, nuclear 
parton distributions are discussed for explaining high-energy nuclear reaction data. Possible nuclear modification 
is explained for valence-quark and antiquark distributions. It is rather difficult to determine gluon distributions in 
nuclei. Next, reversing the topic, we discuss structure functions which could be investigated by neutrino reactions. 
Determination of polarized parton distributions in the nucleon is discussed in polarized neutrino reactions. In 
addition, neutrino reactions should be important for finding nuclear modification of valence-quark distributions 
at small x if structure function ratios / are measured for various nuclei. 



1. Introduction 

Corrections of nuclear effects become impor- 
tant for understanding neutrino cross sections, 
as their measurements become more accurate, in 
order to find the details of neutrino properties 
jl]. In this paper, the nuclear corrections are dis- 
cussed in the high-energy region in terms of par- 
ton distribution functions (PDFs). Unpolarized 
PDFs in the nucleon have been investigated for 
many years, and they are well established from 
very small x to relatively large x. However, it 
is known that parton distributions in a nucleus 
are different from those in the nucleon. Such nu- 
clear modification was first found by the Euro- 
pean Muon Collaboration (EMC) at medium x. 

Now, detailed nuclear modification is known in 
the structure function F2 from small x to large x. 
In the first part of this paper, we discuss possible 
nuclear parton distributions which are extracted 
from deep- inelastic-scattering (DIS) experimental 
data partly with Drell-Yan data j^]. The nu- 
clear PDFs are defined at fixed with a num- 
ber of parameters, which are determined by a 
analysis of the data. There are also other anal- 
yses m Ref. g by including the F^/F^ and 
Drell-Yan data together with dependent data 
ofFf'Vi^f. There are differences between these 
two groups, especially in the x dependent func- 
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tional form of the initial PDFs. From the present 
data, it is rather difficult to determine nuclear 
gluon distributions; however, the dependent 
data could impose some restrictions through 
evolution equations We discuss an analysis 
method and results, from which useful computer 
subroutines are supplied H] for calculating the nu- 
clear PDFs numerically at given x and points 
for a certain nucleus. The obtained nuclear PDFs 
could be used for application to high-energy nu- 
clear reactions such as heavy-ion reactions and 
neutrino- nucleus reactions which are the top- 
ics under discussion at this workshop. 

To the contrary, high-energy neutrino reactions 
should play an important role in determining the 
nuclear PDFs if accurate data are taken for struc- 
ture function ratios to the one for the deuteron. 
However, it is not possible at this stage due to 
lack of accurate deuteron data. In future, there 
are possibilities to have neutrino factories, and we 
discuss possible PDF studies at such factories in 
the last part of this paper. In particular, we ex- 
pect that the valence-quark distributions at small 
X could be investigated. 

Future polarized neutrino measurements 
should be valuable for the studies of nucleon 
spin. It was revealed that quarks carry a small 
fraction of the nucleon spin in contradiction to a 
naive quark model prediction. The neutrino re- 
actions have an advantage of measuring the spin 



1 



content directly. In addition, polarized valence- 
and strange-quark distributions should be inves- 
tigated. 

The nuclear PDFs are explained in Sec.^ 
for understanding nuclear corrections. Then, 
neutrino-nucleon and neutrino-nucleus reactions 
are discussed for better determination of PDFs 
in Sec.^. Our studies are summarized in Sec.|^. 

2. Parton distribution functions in nuclei 

Nuclear medium effects have been investigated 
mainly in the structure function i^2- The ra- 
tio F2 / F2 is about 15% smaller than the one 
for a medium size nucleus at medium x, and it 
tends to increase at large x (~0.9). At small x, 
the ratio F2 / decreases as x becomes smaller, 
which indicates nuclear shadowing phenomenon. 
A different physics mechanism, such as binding, 
Fermi motion, possible subnucleon effect, or mul- 
tiple scattering, contributes to the modification 
in each x region |^. Here, we do not address 
ourselves to such physical mechanisms, and we 
focus on nuclear parton distributions which are 
extracted from deep inelastic experimental data 
and Drell-Yan processes 

Nuclear modification of the PDFs is typically 
less than 20% for medium size nuclei. Because 
they are not much different from those in the 
nucleon, it is technically easier to investigate x 
dependent functional form of the nuclear modifi- 
cation rather than absolute nuclear PDFs. Fur- 
thermore, experimental nuclear structure func- 
tions are often shown by the ratio to the one for 
deuteron such as F^ jF^ ■ Therefore, we define a 
weight function Wi{x, A, Z), which indicates the 
nuclear modification, by 

ff{x,Ql)^w,ix,A,Z)n{x,Ql), (1) 

where A and Z are the mass number and the 
atomic number, /^(xjQq) is the type-i parton 
distribution at Bjorken variable x and Qq in a 
nucleus, and fi{x,Q'^) is the corresponding dis- 
tribution in the nucleon. Because a wide variety 
of date are not available in the nuclear case, fla- 
vor separation of antiquark distributions is not 
possible at this stage. Therefore, as the parton 
type «, we take f^—u^, d^, q^, and g"^ by as- 



suming flavor symmetric antiquark distributions: 
= J-^ = (= q"^). As one of the widely used 
distributions in the nucleon, the leading-order 
(LO) distributions of the MRST parametrization 
are used in our numerical analysis. 

The weight functions are parametrized and 
they are determined by high-energy nuclear re- 
action data. The nuclear modification part is as- 
sumed to have the following functional form: 

ai{A, Z) + bjX + Cix"^ + djx'^ 

where a^, 6^, Ci, di, and (3i are parameters, which 
are determined by a analysis. In the first anal- 
ysis the 1/A^^^ dependence is assumed for 
simplicity by considering volume and surface type 
contributions to the cross section. We will deter- 
mine this factor from the experimental data in 
future analyses. 

The nuclear PDFs are defined at fixed (= 
Qq), whereas the data are taken at various 
points. The initial point is selected as Q§ = 1 
GeV^. In order to make a analysis with the 
data, the parton distributions are evolved to the 
experimental points by the ordinary DGLAP 
evolution equations The LO analysis is made, 
so that the expression of F^ is give by 

where q^ and q'^ are the quark and antiquark 
distributions in the nucleus A, and Cq is the quark 
charge. The total is given by 

/ Tydata r)theo\2 
j ^ j ' 

The notation R is defined by the structure- 
function ratio R = F^/F^ where D denotes 
the deuteron. In Ref. [pi, Drell-Yan cross-section 
ratios are added to the above definition. 

There are still many parameters to be deter- 
mined by the data. The nuclear distributions 
should be constrained by the three conditions: 
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charge, baryon-number, and momentum conser- 
vations. They are expressed by the parton distri- 
butions as 

Charge: Z = J dx ^ {2 - d^) , 
Baryon A = fdx^{u^ + d^) , (5) 



dx J {u^ + d^) , 
Momentum: A= / dx Ax {u{^ + d^ - 



6q^ + g^) 



where the parton distributions are given by those 
per nucleon as obvious from Eqs. (0) and (||). 
Therefore, three parameters can be fixed by these 
conditions. Other parameters are determined so 
as to minimize the total in Eq. (|^). 

The analysis with the weight functions in Eq. 
is called a cubic polynomial type. Another 
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Figure 1. Comparison with the lithimn data. The 
curve indicates the cubic analysis result at 
GeV2. 





Figure 2. Comparison with the calcium data. 



Figure 3. Comparison with the gold data. 



analysis was also done without the diX^ term, and 
it is called a quadratic polynomial type. Because 
quadratic results are similar to the cubic ones [D , 
they are not shown in the following discussions. 
As a result of the cubic type analysis, the Xmin 
per degrees of freedom becomes x^j„/d.o./.=1.82 
for the 309 total data points. Obtained distribu- 
tions are compared with some of the used data 
in Figs. |l|, H, and ||. In these figures, the solid 
curves are the fit results in the cubic analysis at 
GeV^. Because the data are taken in gen- 
eral at different points, the curves cannot be 
directly compared with the data. However, con- 
sidering that the dependence of the ratio R 
is not fairly large, we find that the curves agree 
with the experimental data from small x to large 

X. 

From the analyses, we obtain the weight func- 
tions for the calcium nucleus in Fig. ^. By defini- 
tion, the curves indicate the nuclear modification 
at Q^ = l GeV^. The valence-quark distributions 
show the depletion at a; ~ 0.6 so as to explain 
the EMC effect of F2 at medium x, and they in- 
crease at large x to explain the Fermi-motion part 
of the F2 data. Because the valence distributions 
do not contribute to F2 at small x, their determi- 
nation is difficult. As discussed in the next sec- 
tion, if accurate neutrino data could be taken for 
the ratio F^/F^ in future, the determination be- 
comes possible. The nuclear gluon distributions 
are difficult to be determined in the whole x re- 
gion at least at this stage. Because the antiquark 
distributions dominate the structure function F^ 
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preliminary 



Q- = 50 GcV= 



Figure 4. Obtained weight functions are shown 
for the valence-quark, antiquark, and gluon dis- 
tributions in the calcium nucleus. 



Figure 5. Experimental Drell-Yan cross section 
ratios are compared with the calculated antiquark 
ratio q^'^/q^. 



at small x, their determination is possible in the 
small X region. However, they cannot be fixed in 
the medium and large x regions. 

Drell-Yan data are valuable for determining 
the antiquark distributions, and their analy- 
sis is in progress The Drell-Yan measure- 
ments are shown by the ratios (^^oy/ct'oy ^^'^ 
"'St'/'^ZJY^'^^^- Parton momentum fractions are 
expressed xi and X2 in the proton and a nu- 
cleus, respectively. In the large xp i=xi —X2) re- 
gion, namely in the small X2 region, the Drell-Yan 
cross-section ratio cr^'y/cf^y' is roughly propor- 
tional to the antiquark-distribution ratio q^/q^ . 
A typical result is shown in Fig. ^ together with 
the experimental data. The dashed curve is the 
antiquark ratio q^^/q^ at Q^=50 GeV^. It is 
almost equal to the cross section ratio cTiyy /'^dy 
except for the large X2 region (x2 > 0.2). The the- 
oretical curve is based on our preliminary analysis 
results. 

As shown in the figure, the antiquark distribu- 
tions are restricted in the region x ^ 0.1, where 
there is almost no nuclear modification. There- 
fore, the Drell-Yan data are valuable for finding 
the antiquark distributions in the medium x re- 
gion 01 . We hope to report the analysis results 
with the Drell-Yan data in the near future 

From the analysis results of Ref . |^ , we pro- 
vide subroutines for calculating the parton distri- 
butions at any x and for a given nucleus. The 



codes could be obtained from our web site [g| . Al- 
ternatively, one could use the obtained analytical 
expressions of the weight functions Wi{x, A, Z) for 
calculating the parton distributions. 

There are other analyses by Eskola, Honka- 
nen, Kolhinen, Ruuskanen, and Salgado [Q, and 
their results are rather different from the ones dis- 
cussed in this section. First, for the initial distri- 
butions, namely the weight functions in Eq. (|^), 
they divided the x region into three parts: shad- 
owing, EMC, and Fermi- motion regions. Then, a 
different x dependent functional form is employed 
in each x region for fitting the DIS experimental 
results with the Drell-Yan data. However, as far 
as we are aware, the details of the analysis, es- 
pecially values, have not been reported yet. 
In our first version the Drell-Yan data were 
not included in the analysis partly because we 
intended to find how well the PDFs could be de- 
termined solely by the F^/F^ data. Therefore, 
the antiquark modification increased monotoni- 
cally as x becomes larger in contrast with the 
one in Ref. where the modification is sup- 
pressed in the x ^ 0.1 region. In the same way, 
the gluon distributions at medium x are also dif- 
ferent; however, they could not be determined 
well in any case in the present situation. Obvi- 
ously, these analyses are still premature, so that 
nuclear parametrizations should be investigated 
further by technical refinements with new exper- 
imental data. 
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3. Structure functions in neutrino reac- 
tions 

In the previous section, we discussed the deter- 
mination of the nuclear PDFs by using the elec- 
tron and muon DIS data partly with the Drell- 
Yan data. In this section, we discuss the role of 
high-energy neutrino scattering in the determina- 
tion of the PDFs in the nuclcon and nuclei. 

3.1. Polarized parton distributions 

Neutrino reaction data have been used for ex- 
tracting the unpolarized PDFs in the nucleon 
although actual targets are mainly the iron. 
They are valuable especially for determining the 
strange- and valence-quark distributions. Now, 
the major part of the unpolarized PDFs is well 
known, and our interest tends toward polarized 
PDFs in the nuclcon. Available experimental in- 
formation is limited for the polarized PDFs in the 
sense that the data come mainly from inclusive 
electron and muon DIS and some semi-inclusive 
processes. Furthermore, the small x {x ^ 10^*) 
part is not measured yet, which makes it difhcult 
to determine polarized antiquark distributions. 
This fact affects the determination of the quark 
spin content AS, which has been controversial for 
a decade. Future polarized neutrino reaction data 
should be able to provide important information 
on the polarized distributions and the quark spin 
content. 

In polarized neutrino-proton reactions, struc- 
ture functions g^^ and f/J'^ are expressed by the 
parton distributions M: 



=Ad + As + Au + Ac, 



9i 



~Au + Ac + Ad+ As, 



(6) 



in the leading order (LO). Combining these ex- 
pressions, we find 



dx (g^^ 



97) 



AS. 



(7) 



Although there are higher-order corrections, the 
integral is directly proportional to the quark spin 
content AS of the nucleon in the LO. This is 
an advantage of the neutrino reactions in com- 
parison with the electron and muon scattering, 



where AS cannot be observed directly due to the 
quark charge factor e^. In fact, the obtained spin 
content from the present DIS data ranges from 
AS=0.05 to 0.3. This uncertainty comes from 
the difficulty in determining the polarized anti- 
quark distributions at small x jiot . 

In addition to this advantage, the neutrino 
reactions provide valuable information by the 
parity- violating structure functions gs, g^, and 175 
In the parton model, they are expressed 



{gr + gr)/2x = g7 

= ~{Ad + As- Au- Ac), 
(fff + 5f)/2x = .9f^ 

= -(Am + Ac- Ad- As). (8) 

Combining these structure functions, we have 

97 + 9? = -{^u^ + Ad„) - (As - As) 
-(Ac- Ac), 

I [53'^'+"^ - .93'"+"^ ] - (As + As-) - (Ac -I- Ac). 



(9) 



Because the differences As — As and Ac — Ac are 
expected to be small, the combination g^^ + g^^ 
could be used for finding the polarized valence- 
quark distributions. On the other hand, the pro- 
ton and neutron combination 53 — 33 ^^""""^ is 
valuable for determining the strange and charm 
quark distributions. Furthermore, the strange- 
quark polarization As could be investigated by 
the opposite-sign dimuon events. In this way, 
the neutrino scattering should be able to provide 
valuable information on the polarized parton dis- 
tributions, which cannot be probed by the elec- 
tron scattering and hadron-hadron reactions. 

3.2. Parton distributions in nuclei 

We have discussed the nuclear parton distribu- 
tions in Sec. |[ We found that their determination 
is not straightforward by the electron and muon 
scattering data even together Drell-Yan cross sec- 
tions. Neutrino-nucleus reactions should be able 
to provide valuable information. Actually, there 
are many available neutrino-iron scattering data. 
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However, because there is no accurate neutrino- 
deuteron DIS data, it is almost impossible to in- 
vestigate the nuclear modification at this stage. 
If an intense neutrino beam becomes available 



at a neutrino factory |12|, it becomes possible 
to investigate neutrino reactions with the proton, 
deuteron, and various nuclei. 

It is especially interesting to investigate the 
structure function F3, which docs not exist in 
the electron and muon scattering. A combina- 
tion of neutrino and anti-ncutrino cross-sections 
is expressed in the parton model as 



1 



Fl 



+ + {s ~ s) + (c - c) 



(10) 



The differences between the antiquark distribu- 
tions are considered to be small, so that the struc- 
ture function F3 is appropriate for probing the nu- 
clear modification of the valence-quark distribu- 
tions. Because the ratio / F^ should be equal 
to F^ /F2 in the LO at medium and large x, the 
advantage of the neutrino reactions is to inves- 
tigate the valence-quark modification at small x. 
For example, the analysis results in Sec. || indicate 
very small shadowing at small x as shown by the 
solid curve in Fig. 0. On the other hand, some 
shadowing models produce the modification 
which is similar to the one for the structure func- 
tion i^2, and it is shown by the dotted curve. It 



indicates strong shadowing in the valence-quark 
distributions. Because the EMC effect of F2 at 
medium x should be explained by the valence- 
quark modification, both curves have depletion at 
X 0.5. Due to the baryon number conservation 
in Eq. (||), the modification should be positive 
in the region, x 0.1. However, the valence- 
quark distributions are small in the small x re- 
gion, so that they are not strongly constrained 
by the baryon number. At this stage, it is not 
even obvious whether they show shadowing be- 
havior or antishadowing at small x. It is very im- 
portant to determine the nuclear valence-quark 
distributions at small x in order to find accurate 
antiquark and gluon distributions in nuclei and to 
test nuclear shadowing models. However, there is 
no information on the valence-quark shadowing 
at this stage because the antiquark distributions 
dominate the structure function F2 at small x. 

In addition, if accurate experimental data are 
taken on neutrino- nucleus reactions, it becomes 
possible to investigate nuclear antiquark distri- 
butions by combing different structure functions 
in the similar way to the polarized case of Eq. 

. A future neutrino factory should play an im- 
portant role for the determination of the nuclear 
parton distributions. Furthermore, the details of 
antiquark flavor asymmetry and isospin violation 
could be investigated in the nucleon and nuclei 
by combining different structure functions in the 
neutrino reactions jl^ . 
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Figure 6. Nuclear modification of the valence- 
quark distributions. 



4. Summary 

Nuclear parton distribution functions are de- 
termined by high-energy nuclear reaction data, 
so that they could be used for calculating nu- 
clear corrections in high-energy neutrino reac- 
tions. Nuclear shadowing at small x, the EMC 
effect at medium x, and Fermi motion at large x 
are well explained by the parametrization. How- 
ever, precise determination is still difficult in some 
X regions. For example, it is not straightforward 
to obtain the valence-quark distributions at small 
X and the antiquark distributions at medium x. 
The gluon distributions could not be fixed well 
in the whole x region. The nuclear PDF studies 
are still at the preliminary stage; therefore, future 
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theoretical and experimental efforts are crucial for 
the precise determination. 

On the other hand, if a neutrino factory is 9. 
built in future, the details of the parton distri- 
butions in the nucleon and nuclei could be in- 10. 
vestigated. We showed that the valence-quark 
shadowing should be clarified by the measure- 
ments of structure functions for various nuclei. 11. 
Furthermore, polarized neutrino reactions should 
clarify the issue of quark spin content in the nu- 
cleon by the measurements of the spin dependent 12. 
structure functions. The polarized measurements 
should lead to better determination of the polar- 
ized PDFs. 13. 
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